Finite-Time Euler singularities: A Lagrangian perspective 
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We address the question whether a singularity in a three-dimensional incompressible inviscid 
fluid flow can occur in finite time. Analytical considerations and numerical simulations suggest 
high-symmetry flows being a promising candidate for a finite-time blowup. Utilizing Lagrangian 
and geometric non-blowup criteria, we present numerical evidence against the formation of a finite- 
time singularity for the high-symmetry vortex dodecapole initial condition. We use data obtained 
from high resolution adaptively refined numerical simulations and inject Lagrangian tracer particles 
to monitor geometric properties of vortex line segments. We then verify the assumptions made by 
analytical non-blowup criteria introduced by Deng et. al [Commun. PDE 31 (2006)] connecting 
vortex line geometry (curvature, spreading) to velocity increase to rule out singular behavior. 
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Introduction The incompressible Euler equations in 
three dimensions arc 



du „ 
— + u • Vu 

at 



Vp = , V • u = . 



(1) 



Existence and uniqueness of its solutions for all times 
are unknown. Together with its prominent brother, the 
incompressible Navier-Stokes equations, these equations 
have withstood the minds of mathematicians and physi- 
cists for centuries. While the latter are regarded as "Mil- 
lennium Prize Problem" by the Clay Mathematics In- 
stitute [lj], the ignorance regarding existence of global 
solutions is even larger for the inviscid case: The no- 
tion of weak solutions, which is well established for the 
Navier-Stokes equations since Leray Q, is unknown for 
the three-dimensional Euler equations. 

As a now classical result, the blowup criterion of Beale 
ct al. Q (BKM) connects the existence of solutions for 
the incompressible Euler equations in three dimensions to 
the critical accumulation of vorticity. It has been tried in 
the past to construct explicit initial conditions to obtain 
numerical evidence for or against a finite-time singular- 
ity via BKM, with surprisingly inconsistent results 0, 0] . 
The major reason for this ambiguity is the critical depen- 
dence on extrapolation, which renders the identification 
of singular versus near-singular behavior next to impossi- 
ble by numerical means. The hopes are high that the sit- 
uation is less vague when considering geometric analysis 
of the flow In this letter, we present the application 

of such geometric criteria to numerical data to sharpen 
the distinction between singular and near-singular flow 
evolution. 

The Letter is organized as follows: we first review the 
notion of geometric non-blowup criteria and state the 
considered criteria and their interpretation. We then 
briefly name the computational setup and implementa- 
tion details of our numerical scheme to obtain adaptively 
refined data of up to 8192 3 mesh points. Using La- 
grangian tracers and diagnostics for vortex line geometry, 



we analyze the simulation data to conclude a non-blowup 
for the considered initial conditions. A conclusion and 
outlook summarize the Letter. 

Geometric non-blowup criteria. Historically, non- 
blowup criteria for the incompressible Euler equations 
commonly focus on global features of the flow, such as 
norms of the velocity or the vorticity fields. This comes 
at the disadvantage of neglecting the structures and phys- 
ical mechanisms of the flow evolution. A strategy to 
overcome such shortcomings was established by focus- 
ing more on geometrical properties and flow structures, 
such as vortex tubes or vortex lines. Starting with the 



works of Constantin 
doba and Fcffcrman 



1 01 ] , Constantin et al. |6[ , and Cor- 
J[ j some of these geometric criteria 
(e.g. [8, 11, HI) have reached a phase where they allow 
direct verification of their assumptions with the help of 
numerical simulations. 

Common to geometric criteria is the notion of vortex 
lines, defined as integral curves along the vorticity direc- 
tion field They are transported with the flow, i.e. two 
points x and y on the same vortex line c(s) stay on the 
same vortex line indefinitely. As simple consequence of 
the solenoidality of u> and the BKM theorem, one gets: 

Deng-Hou-Yu theorem 1: Let x(t) be a family of 
points such that for some cq > it holds \u)(pz(t), t)\ > 
coQ(t). Assume that for all t G [0,T) there is another 
point y(t) on the same vortex line as x(i), such that the 
direction of vorticity £(x, t) = w(x, i)/|u;(x, t)\ along the 
vortex line c(s) between x(£) and y(t) is well-defined. If 
we further assume that 



y(t) 



(V-£) (c(s),t)ds 



x(t) 



< c 



for some absolute constant C, and 

i-T 

|w(y(t),t)|dt<oo, 
then there will be no blowup up to time T. 



(2) 



(3) 
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FIG. 1: For the position x(i) of maximum vorticity, choose 
y(t) such that | (V ■ £) (c(a),t)da j = C. For a point- 
wise singularity, x(t) and y(t) must collapse in finite time. 





X{L tl , £l, £2) 



FIG. 2: Lagrangian evolution of a vortex line segment L t in 
the context of theorem 2. For every t%>t\, choose L t2 such 
that it is a subset of X(Lt 1 , ii, £2). 



This criterion can readily be applied to numerical sim- 
ulations. On the other hand, the same theorem may be 
interpreted in a different way to distinguish between a 
point- wise blowup and the blowup of a complete vortex 
line segment, as sketched in Fig. [1] At each instance 
in time, identify the point of maximum vorticity as x(i). 



Now define y(t) such that 
for a constant threshold C 



(c(s),t)ds 



c 



If a singularity occurs, then 
either y(t) approaches x(i) (point-wise blowup), or, if 
the distance between x(t) and y(i) stays finite, the com- 
plete vortex line segment between x(i) and y(t) exhibits 
critical growth. 

Results obtained with this method can be further im- 
proved by considering the Lagrangian evolution of vortex 
line segments L t in time. The geometric equivalent of the 
vortex stretching term is the increase in length for a La- 
grangian vortex line segment. It is possible to quantify 
this stretching and establish a sound connection to the 
vorticity dynamics of the flow. 

Denote with l(t) the length of a vortex line segment 
L t at time t and define with Vti,(f) := ||u>(-, t) \\L°°(L t ) ^ ne 
maximum vorticity on the vortex line segment. Further- 
more, let M(t) := max(||V • £\\ L °°(L t ), ||«IU°°(j t )) be the 
quantity of vortex line convergence V • £ and vortex line 
curvature k, and define A(L t ) := M(t)l(t). It then holds 
(compare Fig. [2]): 

Deng-Hou-Yu theorem 2: Assume there is a fam- 
ily of vortex line segments L t and Tq £ [0,T), such that 
L t2 C X(L tl ,tx,t 2 ) for all T < h < t 2 < T. We 
also assume that Q(t) is monotonically increasing and 
\\u(t)\\L°°(L t ) — c Q^{t) f or some Co > when t is suffi- 
ciently close to T . Furthermore, we assume that 

(i) U&) + U n (t)X(L t ) < (T- t)- A for A e (0, 1) 

(ii) \(Lt) < Co, 

(iii) Kt) ^(T-t) B for some B < 1 - A. 

Then there will be no blowup in the 3D incompressible 
Euler flow up to time T. 



Here, a(t) < b(t) means that there exists a constant 
c£l such that \a(t)\ < c\b(t)\. The velocity compo- 
nents are defined as U^(t) := max Xiy6 L f |(u-^)(x, t) — (u- 
an( i U„(t) := max^ t |u • n|. The proof is given 

in [8(. 

All information about the geometric properties of the 
vortex line segment under consideration is encoded in 
A(-Lt), characterizing the geometric "tameness" of the 
vortex line filament. As depicted in Fig. |3j a vortex 
line segment has a huge X(L t ), if its maximum curvature 
is large, relative to its length (the segment is "kinked" 
instead of "curved"), or if the surrounding vortex lines 
collapse to the considered segment in at least a point (the 
surrounding is "tightening" instead of "parallel"). This 
quantifies the constricted notion of "relatively straight" 
and "smoothly directed" given in Q in a sharper way. 
The process of keeping X(L t ) bounded by some constant 
Co translates to words as the process of "zooming in" 
to the location of maximum vorticity in order to keep 
the considered vortex line segment relatively straight in 
comparison to its length. The assumed accompanying 
collapse in length to keep \(L t ) bounded is then linked 
in its growth rate to the blowup of the velocity compo- 
nents. 

Despite high hopes from an analytical point of view 
that these considerations will shed light on the true na- 
ture of vorticity accumulation, numerical results observ- 




FIG. 3: Characterizing vortex line geometry in terms of 
X(L t ). A slightly curved vortex line with approximately par- 
allel neighboring vortex lines (left) exhibits small \{L±). Vor- 
tex lines with tightening neighboring vortex lines (center) or 
vortex lines with high curvature, in comparison to their length 
(right) have large X(Lt). 
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FIG. 4: Left: Volume plot of the vorticity for the vortex 
dodecapole initial conditions: Six pairs of anti-parallel vortex 
tubes. Right: Isosurface of the vorticity for a single tube at 
a late time. 



ing geometrical properties of Lagrangian vortex filaments 
are scarce. This is primarily due to the fact that Eu- 
lerian quantities such as f2(t) are readily trackable in 
post-processing, while monitoring the Lagrangian evo- 
lution requires additional computational effort. On top 
of that, the geometry of integral curves at an instance in 
time, though in principle computable in post-processing, 
as well as derived quantities such as their convergence 
and curvature, are quite inaccessible in comparison to 
simple Eulerian criteria. 

Numerical Experiment Vorticity-strain coupling is 
the favored mechanism for the formation of a finite-time 
Euler singularity. Since it is well established that this 
process is inherently unstable for turbulent flows, it seems 
natural to search for techniques to artificially keep the 
coupling existent. One such technique is the introduc- 
tion of symmetries to the flow. Early examples such as 
the Taylor-Green vortex [3] or Kerr's initial conditions 
are employing such symmetries. Yet, as pointed out 
by Pelz 14j, for a single vortex tube to exhibit critical 
vorticity-strain coupling, its curvature in the symmetry 
plane has to blow up alongside the axial strain. On the 
other hand, increasing axial strain diminishes the curva- 
ture of the critical vortex line. These counteracting pro- 
cesses constitute an intrinsic resistance of a single vortex 
line to locally "self-stretch" in a critical way. The same 
argument holds for a pair of anti-parallel vortex tubes. A 
popular way to counter this effect is to induce the axial 
strain by neighboring tubes instead of relying on a suffi- 
ciently large kink. This was accomplished, as suggested 
in |l5j . by introducing additional rotational symmetries 
to the flow, arriving at the vortex dodecapole initial con- 
dition pictured in Fig. [4] (left). These initial conditions 
are recognized as a promising candidate for the forma- 
tion of a finite-time singularity. As additional benefit, 
the high symmetry introduces huge savings in computa- 
tional effort and memory requirements. An isosurface of 
the vorticity of one of the vortex tubes at a late time is 
depicted in |4] (right), showing the typical roll- up of the 
vortex sheet in the critical region. 

Resolution is paramount for a reliable statement on 
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FIG. 5: Length of the vortex lines starting at position x of 
maximum vorticity for different constants C — f V • £ds. 



possible singular behavior of the Euler equations. We 
use adaptively refined meshes provided by the frame- 
work racoon [16| to reach effective resolutions of up to 
8192 3 mesh points, scaling close to optimal on up to 
10 5 cores on massively parallel machines. The numeri- 
cal scheme consists of a strong stability preserving third 
order Runge-Kutta [13] time integrator combined with 
a third order shock-capturing CWENO scheme [l8| to 
reduce oscillations in the presence of strong gradients. 
The integrated equation is the vorticity formulation of 
the Euler equations, employing a vector potential formu- 
lation AA = — u> with u = V x A to ensure solenoidality 
of the vorticity vector field lj. The associated Poisson 
equation is solved with a second order parallel and adap- 
tive multigrid algorithm. Interpolation on the coarse-fine 
interfaces is done in us to ensure the highest possible accu- 
racy when applying the aforementioned blowup-criteria. 
Passive tracer particles are injected into the flow for the 
tracking of Lagrangian vortex line segments. The above 
third order Runge-Kutta is also used for the time inte- 
gration of the tracer particles and the space integration 
of vortex lines. Details of the numerical scheme, regard- 
ing its implementation, adaptivity, parallclization, and 
diagnostics will be presented elsewhere. 

The numerical experiment now consists of two parts: 
First, utilizing the method outlined above, we analyze 
the nature of the possible singularity in terms of locality 
(point-wise versus filament). Then, by means of the sec- 
ond theorem, we conclude from the scaling of length and 
velocity components of the critical Lagrangian vortex fil- 
ament, if the observed behavior is singular at all. 

Results For the first part, the constant C of theo- 
rem 1 is chosen in a reasonable way to achieve a length 
of the vortex line segment that fits into the computa- 
tional domain in the beginning of the simulation, but is 
still well resolved at the chosen resolution at later times. 
Hence, the whole vortex line segment is resolved reliably 
throughout the simulation. 
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FIG. 6: Vorticity at the endpoint y of the considered vortex 
line segments. Once satiated, the growth rate is the same for 
all y. 



The results for the vortex dodecapole arc presented in 
figure [5] for different constants C G {0.25,0.5,1,2}. Ini- 
tially, the vortex line segments do not accumulate enough 
V • £, so that the length is bounded by the size of the 
computational domain (x £ [0,7r] 3 ). At some point, de- 
pending on the value of C, the threshold is reached and 
the length of the vortex line segment decreases. Yet, for 
all considered cases of C, the length does not collapse to 
a point, but saturates at early times without approach- 
ing l(t) = 0. This behavior appears to be stable up to 
the latest time of the simulation. The final length of the 
vortex line segments is at least / = 0.3 for the smallest 
case of C (C = 0.25), which is still well resolved with at 
least 200 Ace for the simulation with 4096 3 grid points. 
This result, therefore, is a numerical evidence against a 
point-wise blowup for the vortex dodecapole class of ini- 
tial conditions. This is in concordance with the estimate 
in 0. 

Yet, monitoring the development of u>(y(t),t), as 
shown in figure [6j yields a similar growth rate for the 
accumulation of vorticity at the endpoint y(t) as for 
the beginning of the vortex line segment x(t). This is 
hardly surprising, since by construction a constant value 
for C directly links the growth rates of |u>(x(t),i)| to 
\u>(y(t), t)\. Nevertheless, a numerical verification of this 
analytic equality may be seen as a confirmation that the 
observed growth rate of |w(x(i), t)\ is by no means a nu- 
merical artifact in an isolated small but is repro- 
duced at points far away from the critical region, which 
appear to be well-behaved at first view. The possibly 
critical growth in the perspective of BKM is, thus, con- 
firmed by the global flow. 

Furthermore, since for a large portion of the simu- 
lation the distance l(t) is approximately constant, this 
could be interpreted as an evidence for the existence of a 
non- vanishing vortex line segment that blows up in every 
point. Thus, contradicting the estimation of Q, the pos- 



sibility of a blowup of the vortex dodecapole flow is not 
excluded by theorem 1. The scenario of a collapse to a 
single point, on the other hand, is clearly conflicting the 
numerical evidence. 

For the second part of the numerical experiment, the 
above mentioned theorem 2 is verified numerically to rule 
out a singularity in finite time. The strategy is as follows: 

• Identify the Lagrangian fluid clement a, which will 
contain the maximum of vorticity at the latest time 
of the simulation, tt(t) « \u(X(at,t),t)\. Numer- 
ically this procedure is implemented by carrying 
out a precursory identical simulation with a huge 
number of tracer particles 1 million) randomly 
distributed across the domain. Particles that accu- 
mulate huge amounts of vorticity are selected for 
the production run. 

• In a subsequent computation, at each instance in 
time, start a vortex line integration at X(a, t) along 
the vorticity direction field. Monitor the maximum 
curvature ||/«||i,=o(L t ) and the maximum vortex line 
convergence ||V • £||L^(L t ) during the integration 
and calculate \(t). Stop the integration, as soon 
as \(t) reaches a fixed, arbitrary constant C . This 
defines L t . 

• For this vortex line segment L t , calculate the length 
l(t), and the velocity components U n and U^. From 
the collapse of the length l(t) approximate the ex- 
ponent B. This in turn provides the critical growth 
exponent A for the velocity variables, A CI a =1—B. 

• Compare the increase in U n and to 1/(T— t) AcTit 
to distinguish between critical and sub-critical 
growth of velocity. 

This can be interpreted rather intuitively. By pre- 
scribing an arbitrarily fixed X(t), the vortex line seg- 
ment is kept relatively geometrically uncritical, as the 
length-scale is always adjusted accordingly. This pro- 
cess of "zooming in" just enough to retain the geometric 
"criticalness" prescribes the rate of collapse to a point, 
at least in the direction of the vortex line. All that is left 
to check is whether the velocity growth in the immedi- 
ate surrounding is fast enough to be compatible with a 
finite-time singularity. 

Figure [7] shows the results for the vortex dodecapole 
initial conditions. Pictured is the length of the vortex 
line segment for the tracer that is arriving at a position 
of very huge vorticity at late stages of the simulation. 
The subplot depicts the long-term behavior of the par- 
ticle entering the critical region, while the final stage of 
length decrease is magnified. The decrease in length does 
not agree with a collapse in final time, but instead the 
shrinkage of the segment decelerates clearly in time. This 
contradicts a scaling in time proportional to (T — t) B for 
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FIG. 7: Evolution of the length l(t) of the critical vortex fil- 
ament L t for different Lagrangian fluid elements. The length 
does not decrease as (T — t) B for any B < 1, which would 
be faster than linear. The Lagrangian collapse of the vortex 
segment is decelerating instead. 




FIG. 8: Evolution of the quantities U n and U% in time. Both 
quantities do not exhibit a critical growth of 1/(T — t). 



any < B < 1, which would be faster than (or, in the 
limiting case, equal to) linear. 

It could furthermore be argued that the limit B — > 
is hard to exclude, since the drop in length would be 
virtually instantaneous in time, with a close to constant 
scaling before. In this limit, the quantities U n and U% 
would have to grow roughly as 1/(T — t) to still allow 
formation of a finite-time singularity. Figure [5] shows 
the observed behavior of U n and in time. Both show 
no signs of critical accumulation, in particular not like 
1/(T — t) in time. Thus, the assumptions of theorem 
2 are well met. These results therefore pose a strong 
evidence against a finite-time singularity for the class of 
vortex dodecapole initial conditions. 



Conclusions and Outlook In this Letter we studied 
the question whether a singularity in a three-dimensional 
incompressible inviscid fluid flow can occur in finite time. 
Using massively parallel high-resolution adaptive simula- 
tions and applying Lagrangian and geometric diagnostics 
to the flow evolution, we are able to probe the behavior of 
the critical vortex filament. Our findings pose a strong 
numerical evidence against a finite-time singularity for 
the vortex dodecapole initial conditions, based on ana- 
lytical criteria connecting velocity scaling to vortex line 
segment geometry In principle, the presented method 
could easily be applied to different classes of initial con- 
ditions. 
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